Missense mutations in the membrane domain of PRRT2 affect its interaction with Nav1.2 voltage-gated sodium channels.

PRRT2 is a neuronal protein that controls neuronal excitability and network stability by modulating voltage-gated Na+ channel (Nav). PRRT2 pathogenic variants cause pleiotropic syndromes including epilepsy, paroxysmal kinesigenic dyskinesia and episodic ataxia attributable to loss-of-function pathogenetic mechanism. Based on the evidence that the transmembrane domain of PRRT2 interacts with Nav1.2/1.6, we focused on eight missense mutations located within the domain that show expression and membrane localization similar to the wild-type protein. Molecular dynamics simulations showed that the mutants do not alter the structural stability of the PRRT2 membrane domain and preserve its conformation. Using affinity assays, we found that the A320V and V286M mutants displayed respectively decreased and increased binding to Nav1.2. Accordingly, surface biotinylation showed an increased Nav1.2 surface exposure induced by the A320V mutant. Electrophysiological analysis confirmed the lack of modulation of Nav1.2 biophysical properties by the A320V mutant with a loss-of-function phenotype, while the V286M mutant displayed a gain-of-function with respect to wild-type PRRT2 with a more pronounced left-shift of the inactivation kinetics and delayed recovery from inactivation. The data confirm the key role played by the PRRT2-Nav interaction in the pathogenesis of the PRRT2-linked disorders and suggest an involvement of the A320 and V286 residues in the interaction site. Given the similar clinical phenotype caused by the two mutations, we speculate that circuit instability and paroxysmal manifestations may arise when PRRT2 function is outside the physiological range.

The impact of pathogenic and artificial mutations on Claudin-5 selectivity from molecular dynamics simulations.

Tight-junctions (TJs) are multi-protein complexes between adjacent endothelial or epithelial cells. In the blood-brain-barrier (BBB), they seal the paracellular space and the Claudin-5 (Cldn5) protein forms their backbone. Despite the fundamental role in brain homeostasis, little is known on Cldn5-based TJ assemblies. Different structural models were suggested, with Cldn5 protomers generating paracellular pores that restrict the passage of ions and small molecules. Recently, the first Cldn5 pathogenic mutation, G60R, was identified and shown to induce Cl--selective channels and Na+ barriers in BBB TJs, providing an excellent opportunity to validate the structural models. Here, we used molecular dynamics to study the permeation of ions and water through two distinct G60R-Cldn5 paracellular architectures. Only the so-called Pore I reproduces the functional modification observed in experiments, displaying a free energy (FE) minimum for Cl- and a barrier for Na+ consistent with anionic selectivity. We also studied the artificial Q57D and Q63D mutations in the constriction region, Q57 being conserved in Cldns except for cation permeable homologs. In both cases, we obtain FE profiles consistent with facilitated passage of cations. Our calculations provide the first in-silico description of a Cldn5 pathogenic mutation, further assessing the TJ Pore I model and yielding new insight on BBB's paracellular selectivity.

Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels.

The recent determination of cryo-EM structures of voltage-gated sodium (Nav) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Nav pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Nav channels based on homology modeling of the cryo-EM structure of the human Nav1.4 channel and, in addition, on the recently resolved configuration for Nav1.2. In particular, single Na+ permeation events during standard MD runs suggest that the ion resides in the inner part of the Nav selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na+ translocation through the SF of the homology-based Nav1.2 model and the cryo-EM Nav1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na+ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Nav cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Nav channels.

The intramembrane COOH-terminal domain of PRRT2 regulates voltage-dependent Na+ channels.

Proline-rich transmembrane protein 2 (PRRT2) is the single causative gene for pleiotropic paroxysmal syndromes, including epilepsy, kinesigenic dyskinesia, episodic ataxia, and migraine. PRRT2 is a neuron-specific type-2 membrane protein with a COOH-terminal intramembrane domain and a long proline-rich NH2-terminal cytoplasmic region. A large array of experimental data indicates that PRRT2 is a neuron stability gene that negatively controls intrinsic excitability by regulating surface membrane localization and biophysical properties of voltage-dependent Na+ channels Nav1.2 and Nav1.6, but not Nav1.1. To further investigate the regulatory role of PRRT2, we studied the structural features of this membrane protein with molecular dynamics simulations, and its structure-function relationships with Nav1.2 channels by biochemical and electrophysiological techniques. We found that the intramembrane COOH-terminal region maintains a stable conformation over time, with the first transmembrane domain forming a helix-loop-helix motif within the bilayer. The unstructured NH2-terminal cytoplasmic region bound to the Nav1.2 better than the isolated COOH-terminal intramembrane domain, mimicking full-length PRRT2, while the COOH-terminal intramembrane domain was able to modulate Na+ current and channel biophysical properties, still maintaining the striking specificity for Nav1.2 versus Nav1.1. channels. The results identify PRRT2 as a dual-domain protein in which the NH2-terminal cytoplasmic region acts as a binding antenna for Na+ channels, while the COOH-terminal membrane domain regulates channel exposure on the membrane and its biophysical properties.

Multiscale modelling of claudin-based assemblies: A magnifying glass for novel structures of biological interfaces.

Claudins (Cldns) define a family of transmembrane proteins that are the major determinants of the tight junction integrity and tissue selectivity. They promote the formation of either barriers or ion-selective channels at the interface between two facing cells, across the paracellular space. Multiple Cldn subunits form complexes that include cis- (intracellular) interactions along the membrane of a single cell and trans- (intercellular) interactions across adjacent cells. The first description of Cldn assemblies was provided by electron microscopy, while electrophysiology, mutagenesis and cell biology experiments addressed the functional role of different Cldn homologs. However, the investigation of the molecular details of Cldn subunits and complexes are hampered by the lack of experimental native structures, currently limited to Cldn15. The recent implementation of computer-based techniques greatly contributed to the elucidation of Cldn properties. Molecular dynamics simulations and docking calculations were extensively used to refine the first Cldn multimeric model postulated from the crystal structure of Cldn15, and contributed to the introduction of a novel, alternative, arrangement. While both these multimeric assemblies were found to account for the physiological properties of some family members, they gave conflicting results for others. In this review, we illustrate the major findings on Cldn-based systems that were achieved by using state-of-the-art computational methodologies. The information provided by these results could be useful to improve the characterization of the Cldn properties and help the design of new efficient strategies to control the paracellular transport of drugs or other molecules.

Computational Assessment of Different Structural Models for Claudin-5 Complexes in Blood-Brain Barrier Tight Junctions.

The blood-brain barrier (BBB) strictly regulates the exchange of ions and molecules between the blood and the central nervous system. Tight junctions (TJs) are multimeric structures that control the transport through the paracellular spaces between the adjacent brain endothelial cells of the BBB. Claudin-5 (Cldn5) proteins are essential for TJ formation and assemble into multiprotein complexes via cis-interactions within the same cell membrane and trans-interactions across two contiguous cells. Despite the relevant biological function of Cldn5 proteins and their role as targets of brain drug delivery strategies, the molecular details of their assembly within TJs are still unclear. Two different structural models have been recently introduced, in which Cldn5 dimers belonging to opposite cells join to generate paracellular pores. However, a comparison of these models in terms of ionic transport features is still lacking. In this work, we used molecular dynamics simulations and free energy (FE) calculations to assess the two Cldn5 pore models and investigate the thermodynamic properties of water and physiological ions permeating through them. Despite different FE profiles, both structures present single/multiple FE barriers to ionic permeation, while being permissive to water flux. These results reveal that both models are compatible with the physiological role of Cldn5 TJ strands. By identifying the protein-protein surface at the core of TJ Cldn5 assemblies, our computational investigation provides a basis for the rational design of synthetic peptides and other molecules capable of opening paracellular pores in the BBB.

Computational study of ion permeation through claudin-4 paracellular channels.

Claudins (Cldns) form a large family of protein homologs that are essential for the assembly of paracellular tight junctions (TJs), where they form channels or barriers with tissue-specific selectivity for permeants. In contrast to several family members whose physiological role has been identified, the function of claudin 4 (Cldn4) remains elusive, despite experimental evidence suggesting that it can form anion-selective TJ channels in the renal epithelium. Computational approaches have recently been employed to elucidate the molecular basis of Cldns' function, and hence could help in clarifying the role of Cldn4. In this work, we use structural modeling and all-atom molecular dynamics simulations to transfer two previously introduced structural models of Cldn-based paracellular complexes to Cldn4 to reproduce a paracellular anion channel. Free energy calculations for ionic transport through the pores allow us to establish the thermodynamic properties driving the ion-selectivity of the structures. While one model shows a cavity permeable to chloride and repulsive to cations, the other forms barrier to the passage of all the major physiological ions. Furthermore, our results confirm the charge selectivity role of the residue Lys65 in the first extracellular loop of the protein, rationalizing Cldn4 control of paracellular permeability.

Structural Mechanism of ω-Currents in a Mutated Kv7.2 Voltage Sensor Domain from Molecular Dynamics Simulations.

Activation of voltage-gated ion channels is regulated by conformational changes of the voltage sensor domains (VSDs), four water- and ion-impermeable modules peripheral to the central, permeable pore domain. Anomalous currents, defined as ω-currents, have been recorded in response to mutations of residues on the VSD S4 helix and associated with ion fluxes through the VSDs. In humans, gene defects in the potassium channel Kv7.2 result in a broad range of epileptic disorders, from benign neonatal seizures to severe epileptic encephalopathies. Experimental evidence suggests that the R207Q mutation in S4, associated with peripheral nerve hyperexcitability, induces ω-currents at depolarized potentials, but the fine structural details are still elusive. In this work, we use atom-detailed molecular dynamics simulations and a refined model structure of the Kv7.2 VSD in the active conformation in a membrane/water environment to study the effect of R207Q and four additional mutations of proven clinical importance. Our results demonstrate that the R207Q mutant shows the most pronounced increase of hydration in the internal VSD cavity, a feature favoring the occurrence of ω-currents. Free energy and kinetics calculations of sodium permeation through the native and mutated VSD indicate as more favorable the formation of a cationic current in the latter. Overall, our simulations establish a mechanistic linkage between genetic variations and their physiological outcome, by providing a computational description that includes both thermodynamic and kinetic features of ion permeation associated with ω-currents.

Genotype-phenotype correlations in patients with de novo KCNQ2 pathogenic variants.

OBJECTIVE: Early identification of de novo KCNQ2 variants in patients with epilepsy raises prognostic issues toward optimal management. We analyzed the clinical and genetic information from a cohort of patients with de novo KCNQ2 pathogenic variants to dissect genotype-phenotype correlations. METHODS: Patients with de novo KCNQ2 pathogenic variants were identified from Italy, Denmark, and Belgium. Atomic resolution Kv7.2 structures were also generated using homology modeling to map the variants. RESULTS: We included 34 patients with a mean age of 4.7 years. Median seizure onset was 2 days, mainly with focal seizures with autonomic signs. Twenty-two patients (65%) were seizure free at the mean age of 1.2 years. More than half of the patients (17/32) displayed severe/profound intellectual disability; however, 4 (13%) of them had a normal cognitive outcome.A total of 28 de novo pathogenic variants were identified, most missense (25/28), and clustered in conserved regions of the protein; 6 variants recurred, and 7 were novel. We did not identify a relationship between variant position and seizure offset or cognitive outcome in patients harboring missense variants. Besides, recurrent variants were associated with overlapping epilepsy features but also variable evolution regarding the intellectual outcome. CONCLUSIONS: We highlight the complexity of variant interpretation to assess the impact of a class of de novo KCNQ2 mutations. Genetic modifiers could be implicated, but the study paradigms to successfully address the impact of each single mutation need to be developed.

Molecular Dynamics Simulations of Ion Selectivity in a Claudin-15 Paracellular Channel.

Claudins are tissue-specific transmembrane proteins able to form junctions between two cells and regulate the flow of physiological solutes parallel to the cell walls, that is, the paracellular transport. Claudin-15 is highly expressed in the intestine where it forms efficient Na+ channels and Cl- barriers. However, the molecular details of these biological complexes are still unclear. Here, the permeation process of Na+, K+, and Cl- ions inside a refined structural model of a claudin-15 paracellular channel is investigated using all-atom molecular dynamics simulations in a double-bilayer and explicit solvent. One-dimensional potential of mean force (PMF) profiles, calculated using umbrella sampling (US) simulations, show that the channel allows the passage of the two physiological cations while excluding chloride. These features are generated by the action of several acidic residues, in particular the ring of D55 residues which is located at the narrowest region of the pore, in correspondence with the energy minimum for cations and the peak for chloride. We also used the Voronoi-tessellated milestoning method to obtain additional PMF profiles and the permeation timescale of the three ions. The milestoning PMFs agree well with those obtained by US, and the rate calculation reveals that the passage of chloride is almost 30 times slower than that of sodium. Our results are consistent with the known ability of claudin-15 to regulate tight junction selectivity and with the experimentally determined role of the acidic residues. This further validates our structural model and provides insights into the atomistic details of ion transport in paracellular channels that could be shared by other claudin-based architectures.

Interfacing Graphene-Based Materials With Neural Cells.

The scientific community has witnessed an exponential increase in the applications of graphene and graphene-based materials in a wide range of fields, from engineering to electronics to biotechnologies and biomedical applications. For what concerns neuroscience, the interest raised by these materials is two-fold. On one side, nanosheets made of graphene or graphene derivatives (graphene oxide, or its reduced form) can be used as carriers for drug delivery. Here, an important aspect is to evaluate their toxicity, which strongly depends on flake composition, chemical functionalization and dimensions. On the other side, graphene can be exploited as a substrate for tissue engineering. In this case, conductivity is probably the most relevant amongst the various properties of the different graphene materials, as it may allow to instruct and interrogate neural networks, as well as to drive neural growth and differentiation, which holds a great potential in regenerative medicine. In this review, we try to give a comprehensive view of the accomplishments and new challenges of the field, as well as which in our view are the most exciting directions to take in the immediate future. These include the need to engineer multifunctional nanoparticles (NPs) able to cross the blood-brain-barrier to reach neural cells, and to achieve on-demand delivery of specific drugs. We describe the state-of-the-art in the use of graphene materials to engineer three-dimensional scaffolds to drive neuronal growth and regeneration in vivo, and the possibility of using graphene as a component of hybrid composites/multi-layer organic electronics devices. Last but not least, we address the need of an accurate theoretical modeling of the interface between graphene and biological material, by modeling the interaction of graphene with proteins and cell membranes at the nanoscale, and describing the physical mechanism(s) of charge transfer by which the various graphene materials can influence the excitability and physiology of neural cells.

Correction: A refined model of claudin-15 tight junction paracellular architecture by molecular dynamics simulations.

[This corrects the article DOI: 10.1371/journal.pone.0184190.].

A refined model of claudin-15 tight junction paracellular architecture by molecular dynamics simulations.

Tight-junctions between epithelial cells of biological barriers are specialized molecular structures that regulate the flux of solutes across the barrier, parallel to cell walls. The tight-junction backbone is made of strands of transmembrane proteins from the claudin family, but the molecular mechanism of its function is still not completely understood. Recently, the crystal structure of a mammalian claudin-15 was reported, displaying for the first time the detailed features of transmembrane and extracellular domains. Successively, a structural model of claudin-15-based paracellular channels has been proposed, suggesting a putative assembly that illustrates how claudins associate in the same cell (via cis interactions) and across adjacent cells (via trans interactions). Although very promising, the model offers only a static conformation, with residues missing in the most important extracellular regions and potential steric clashes. Here we present detailed atomic models of paracellular single and double pore architectures, obtained from the putative assembly and refined via structural modeling and all-atom molecular dynamics simulations in double membrane bilayer and water environment. Our results show an overall stable configuration of the complex with a fluctuating pore size. Extracellular residue loops in trans interaction are able to form stable contacts and regulate the size of the pore, which displays a stationary radius of 2.5-3.0 Å at the narrowest region. The side-by-side interactions of the cis configuration are preserved via stable hydrogen bonds, already predicted by cysteine crosslinking experiments. Overall, this work introduces an improved version of the claudin-15-based paracellular channel model that strengthens its validity and that can be used in further computational studies to understand the structural features of tight-junctions regulation.